Low dimensional behavior in three-dimensional coupled map lattices 
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The analysis of one-, two-, and three-dimensionaf coupfed map lattices is fiere devefoped under 
a statistical and dynamical perspective. We show that the three-dimensional CML exhibits fow 
dimensionaf behavior with fong range correfation and the power spectrum foffows 1/f noise. This 
approach feads to an integrated understanding of the most important properties of these universal 
models of spatiotemporal chaos. We perform a complete time series analysis of the model and 
investigate the dependence of the signal properties by change of dimension. 
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I. INTRODUCTION 

The widespread interest in coupled map lattices (CML) 
in research and applications is mainly due to their role 
as a family of systems possessing universal behavior 
with explicit form of the local interactions. Applica- 
tions of CML are found in several areas such as dy- 
namics, turbulence, phase transitions, geophysics, op- 
tics, genetics, human information processes, and so forth 
[U, BM i, II i, 0, i, i, El m ■ It is worth to mention that 
Kaneko has studied spatiotemporal pattern dynamics in 
both ID and 2D CML In particular, Kaneko has 

shown the existence of similar patterns of different vari- 
eties in ID and 2D CML in the weak coupling regime. 
However, in the strong coupling regime, the dynamics 
in ID and 2D lattices differ [13] ■ These coupled map 
lattices also exhibit the size instability wherein the dy- 
namics is completely dependent on system size [H, [lj| . 
For example, ID CML shows stable synchronous state at 
small number of lattices while it transits to spatiotempo- 
ral chaos at large system size. 

Although in many instances the focus is on one- and 
two- dimensional lattices, in this paper we will obtain a 
wide range of properties that characterize the time series 
of ID, 2D and 3D coupled map lattices. In this sense a 
comparison can be achieved on the effect of the behavior 
of the system under change of dimension. Due to the 
complexity of this endeavor we will concentrate on the 
logistic maps but the numerical work can be extended, 
using more powerful computers, to Lorenz-type and other 
models. 

As mentioned above, the dynamical dependence on the 
control parameters, particularly in ID and 2D coupled 
map lattices is well known [12, [H [l^ and we fix val- 
ues that allow interesting chaotic regime. In the present 
study, we deal specifically with the time series generated 
by the evolution and explore their statistical, dynamic 
and geometric properties; we also obtain measures of di- 
mension, and predictability when the dimension of the 
lattice is increased. The robustness of our conclusions 



was tested by taking several random points on the lat- 
tice, for each dimension, and analyzing the corresponding 
time series. In spite of the fact that the properties of the 
maps are known, subtle interrelation of the many modes 
of evolution will occur that are not present in ID lattices. 

Several systems of practical interest exhibit 1/f scal- 
ing behavior (TtI . [isj and in this paper we implement de- 
trended fluctuation analysis (DFA) to analyze such prop- 
erty [H, [13, [m . We use statistical and other characteri- 
zations of the CML in order establish a set of properties 
that appear in the time series. Measures of dimension 
and other diagnostics will be used to detect differences 
of behavior of these models. 

In this contribution we aim to provide several argu- 
ments that serve to compare and clarify the relationship 
between the time series of ID, 2D and 3D CMLs. A new 
feature that emerges in 3D models, not present in lower 
dimensions, occurs in the plot of the root mean square 
of the fluctuation. Another is the increase of persistence 
when the dimension is increased. 

The paper is organized as follows. In Section [TT] we 
describe how the model and its numerical implementa- 
tion is built. In Section IIIII we perform the statistical 
analysis, which include DFA and power spectral density. 
The dynamics part of the paper, which includes local di- 
mension analysis, embedding dimension, stationarity and 
predictability, is discussed in Section ITVl In the last sec- 
tion we present a summary of results and conclusions. 



II. DESCRIPTION OF THE MODEL 



Let us first consider an one dimensional coupled map 
Jattices with nearest neighbor coupling defined as [H, [Tj, 



, = {l~e)fixl) + ^[f{xr) + fixl^')], (1) 



where i — 1,2, ... ,N are the lattice sites and N is the 
number of lattice points. Each of the lattice points in 
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Eq. Il]) is represented by the logistic map, 

f{x)^tix{l-x), xe{0,l), AiG(0,4). (2) 

In the two dimensional lattices, each map is coupled to 
four of its nearest neighbors and given by the following 
form: 



{i-e)f{x^^) + '-\f{xr'^^)+f{xr^'-^) 



(3) 



One can also think of a coupled map lattices in three 
dimensions where each of the individual map is coupled 
to six of its neighbors and represented by the following 
set of coupled equations. 

= (1 - e)f{x'^^^'') + 1 \f{xr'^''') + f{x^'''^') 



(4) 



For the present study we fix the parameters as = 4, e = 
0.4 and the size as = 50. The time evolution generated 
by these maps is shown in Fig. [T] The distinguishing 
feature of the 3D case. Fig. [IJc) is the variation of the 
amplitude of the motion resulting in the clustering of 
variances. 
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III. STATISTICAL ANALYSIS 

A. Detrended fluctuation analysis 

Consider a time series a;(z), i = 1, 2, . . . , TV and define 
an integrated the time series y(i) as. 



where x is the mean given by 

N 



(5) 



(6) 



Then the integrated time series y(i) is divided into I boxes 
of equal size, say, n where I = Int {N/n). The data in 
each box is fitted with a linear regression function j/fit(*)i 
the local trend. Now detrend the integrated time series 
by subtracting the local trend 2/fjt(i) in each box. Then 
the root mean square (rms) fluctuation is calculated as 



F{n) = 



/-I ( [k+l)n 

yE E [^w 

A:=0 yi=kn+l 



(7) 



By repeating the above calculation for different box sizes 
the relationship between the rms fluctuation F{n) and 



FIG. 1: (Color online) Time series of a randomly chosen lat- 
tice point for (a) one-, (b) two- and (c) three dimensional 
CML 



the size n is established. In general, there follows a 
power-law relation between F{n) and n, and F{n) scales 
as F{n) ~ n", where a is called scaling exponent or DFA 
exponent or sometimes referred as correlation exponent. 
Depending on the value of a one can classify the time 
series whether it is correlated or anti-correlated. For an 
uncorrelated time series a = 0.5, for positive correlations 
a > 0.5, and a < 0.5 for anti-correlated time series. A 
detailed analysis on the different trends of artificial data 
sets are found in Ref. [2l|. 

When the autocorrelation C{n) separated by n time 
units, decays as 



C(7 



(8) 



with < 7 < 1, the fluctuation exponent obeys the fol- 
lowing relation 
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(9) 



Note that 7 > 1 for uncorrelated data. Further the DFA 
exponent is related to the 1//^ spectral slope as 



f3 = 2a -1, 



(10) 
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and a = 1 for a 1// noise spectrum and a = 1.5 for 
Brownian motion (random walk) where (3 = 2. 

We calculate the exponent a for one-, two- and three- 
dimensional coupled map lattices. In all the three cases 
F{n) scales as a power law. The exponent a for one- 
dimensional CML is found to be 0.54 which is close to 
0.5. This essentially indicates that the time series is un- 
correlated. In the case of two-dimensional CML the ex- 
ponent a ~ 0.62 showing a positive correlation in the 
time series [see Fig. [5] . 
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FIG. 2: (Color online) The variation of F{n) vs. n for the 
one-, two- and three-dimensional CML 

The time series from three-dimensional CML exhibits 
an entirely different behavior. In Fig. [2] we show the 
variation of the fluctuation exponent, F{n), as a function 
of n. It is easy to see from Fig. [5] that there exists two 
regimes with a crossover, which occurs around n ~ 10^. 
Before the crossover there is no apparent scaling present 
in the time series. However, for small n two dominant 
clusterings present in the data. 

For n > 10^, F{n) scales as n" where a « 1 indicat- 
ing a strong positive correlation in the time series. One 
should note that when a = 1, the auto-correlation expo- 
nent 7 = [see eq. ©]. It also indicates that the time 
series follows 1/ noise [see eq. below]. 



B. Power spectral density 

Spectral density is often used to characterize many 
noisy time series. In the present study we use the power 
spectral density (PSD) analysis to find the signature of 
1/f noise as predicted by the DFA analysis. The power 
spectral density for a time series x is defined as 



S{f) - |^(x)|^ 



(11) 



where T{'x.) denotes fast Fourier transform (discrete 
Fourier transform) of the time series x. We are searching 
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FIG. 3: (Color online) Plot of the power spectral density as 
a function of frequency showing 1//'' noise spectrum for the 
time series obtained from 3-dimensional coupled map lattices 



for a relationship of the form 



sif) - Jp 



(12) 



Fig. [3] shows the plot of the PSD as a function of fre- 
quency. Here we take the average of S{f) over a range 
of frequencies A/. From Fig. [3] one can see that, the 
PSD falls as almost l//'^ in the case of three-dimensional 
CML, with /? ~ 1.18. This is consistent with the value of 
as in Fig. [Hand eq. pIT]) . 

Feigenbaum ^] has shown an example of determinis- 
tic dynamical systems with 1/f noise behavior. It has 
also been confirmed theoretically as a real 1// spectrum 
density using simple quadratic map f23l|. It is a known 
fact that there is no generally accepted explanation for 
this phenomenon, in spite of its universality. However, 
in the context of dynamical systems, especially in chaotic 
systems, the intermittency is one of the sources for 1/f 
noise, and as a consequence these systems exhibit long 
range time and spatial correlations f24!| . 

Earlier, we found that the fluctuation exponent for ID 
CML is almost close to 0.5, similar to that observed in 
white noise. For 2D CML the exponent increased slightly. 
However, there is a dramatic change in the exponent for 
3D CML with awl [see Fig. [1]. This indicates the 
signature of 1// noise. The power spectral density of the 
corresponding time series closely follows 1/ f^ decay [see 
Fig. [3] . The structure of correlation on data, according 
Fig. [21 clearly affected the power spectrum producing 
deviation around the 1/f curve. 



C. Structure function analysis 

Structure function (SF) constitute a technique used to 
identify stationarity during short time evolution. It has 
been successfully employed in study in turbulent flows 
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|2a | but it has a much wider applications. The SF is 
defined in the following way: 



S,ir)^{\YiU+r)-YiU)n 



(13) 



where Y{t) denotes some stochastic process, or the com- 
ponent of some dynamical system, i denotes the i-th data 
point, (• • • ) is the ensemble average and g is a real num- 
ber. When q = 2 the SF is just the correlation function. 
liY{t) is scale invariant and self-similar (self-afhne) over 
some time lag interval {ti < t < T2), then the g-th order 
structure function is expected to be fitted by 



(14) 



Note that Cq can be a function of r whose variation is 
varies slower than any power of t, and ((q) is the expo- 
nent of the structure function. 

For stationary time series the exponent of the struc- 
ture function is zero ({q) = 0, as a consequence of the 
translational invariance of any statistics defined from the 
data. Fig. 0] shows the SF for q — 2. For small time lags 
we find two branches meaning that a breakdown of sta- 
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the BV dimension analysis for the 2D and 3D CML. In 
the case of 3D CML we consider a two dimensional slice 
in order to compute the local dimension. The calculation 
of BV dimension is formulated as follows. 

Consider a 2D spatially distributed system whose state 
at a given time ti is defined over a collection of points 
(«, j). Here we take the M —1 nearest neighbors for each 
point (i,j) in a square lattice with M = 25. Logistic 
maps are one variable dynamical systems and in order to 
specify the corresponding state at a point including its 
neighbors we need an M dimensional state vector called 
bred vector. Now generate k distinct perturbations of the 
state starting at to < ti obtaining k local bred vectors. 
The k X k covariance matrix of the system is just C = 
6*^6, where B is the Af x k matrix of local bred vectors 
each normalized to unity. 

We order the eigenvalues of the covariance matrix as 
Ai > A2 > . . . > Afc and define the singular values of B as 
ci — n/A/", I — 1,2, ... k. The eigenvalues A; represent the 
amount of variance in the set of the k unit bred vectors. 
Then the bred-vector dimension can be defined as [261 



i^ij {(^1,0-2, • ■ ■ , cTfc) 



(15) 



As each of the k bred vectors is normalized to unity, ip 
assumes values in the interval (0, k). In the present study, 
we fix the value of the number of bred vectors as fc = 5. 

We compute the local dimensions for the 2D and 3D 
CML for different set of initial conditions. In the case of 
2D CML the local dimensions are found to lie between a 
minimum of 1.03 and a maximum of 2.4, while for the 3D 
CML it takes a minimum value of ~ 1.04 and a maximum 
value of ~ 2.1. Since this kind of dimension is related to 
predictability [13, Ell, we conclude that the 3D CML is 
more predictable than the 2D case. 



FIG. 4: (Color online) Plot of the structure function Sq{T) 
as a function of delay r for the time series obtained from 
3-dimensionaI coupled map lattices 

tionarity in this region. A similar effect has been found 
in the study of the DFA plot in Fig. [2l For higher time 
lags the fiat SF plot is indicative of stationarity in this 
region. 

IV. DYNAMICAL ANALYSIS 
A. Local dimension 

An interesting statistics used to characterize spa- 
tiotemporal chaotic systems is given by bred vector (BV) 
dimension [2^, [l^jllsl- One can well analyse the local di- 
mension of a two-dimensional spatiotemporal chaotic sys- 
tem using BV dimension. In the present study we apply 



B. Embedding dimension 

The method developed in Ref. [2^ was implemented 
to obtain the embedding dimension of the CMLs. In our 
simulations this procedure was more robust and objec- 
tive than the ambiguous results based on the correlation 
integral. The idea is to start with some initial embedding 
dimension d and a vector rj from the set of state vectors. 
A simple prediction algorithm, similar to the one used in 
the cross prediction test of nonstationarity (see below), 
consists of finding the nearest neighbor r^^j^ of rj and 
then to evolve it one time step ahead obtaining rQ,(j)^i. 
The error in this prediction is || rj^i — r^Q).).]^ || and the 
error over all the set of state vectors is just the sum, E{d), 
of each individual error. Then the embedding dimension 
is the value of d that minimizes the function E{d). This 
is a simple and clever way of determining the embedding 
dimension that does not require large data sets. [See, for 
example, Ref. [soj an estimate of the size of a time series 
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that is necessary to give meaningful results when com- 
puting dimensions using the correlation integral. Also in 
Ref. [311 an even more strict estimate is provided]. In 
Fig. [5] we show the dimension of the CMLs. We find that 
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FIG. 5: (Color online) Estimation of optimal embedding di- 
mension for the time series obtained from one-, two- and 
three-dimensional CMLs 

the ID CML has the lowest dimension followed closely 
by the 2D and 3D lattices. All of them have dimensions 
compatible with chaotic systems. 






FIG. 6: (Color online) Recurrence plot for the (a) nonstation- 
ary baker's map H, (b) ID CML, (c) 2D CML and (d) 3D 
CML. 



C. Stationarity and recurrence analysis 

One of the most important tests of dynamical systems 
is related to stationarity. We implement three differ- 
ent algorithms since there is no single criterion for these 
tests in nonlinear dynamics. We also implement a recent 
method [s^j using recurrence analysis to find whether 
these systems possess a positive maximal Lyapunov ex- 
ponent. 

The recurrence plot was introduced in Ref. [33| and 
it provides insights about many properties of dynamical 
systems in a two dimensional graphic display. The idea 
is to plot the times when the trajectory visits a given 
region in phase space. Deterministic nonlinear dynami- 
cal systems exhibit the fundamental property that states 
tend to return to the neighborhoods of regions already 
visited by the trajectory. 

Consider the following function 

i?y = 6(6- II - r, II) (16) 

where the state vectors r^, rj are considered at discrete 
times with I < i, j < N , N the total number of states. 
Here Q is the Heaviside function and e a positive number 
whose value is discussed below. A graph is then com- 
posed by associating a black dot in the plane with co- 
ordinates whenever the argument of the function is 
positive. We examine the well known baker's map 34] for 
comparison and also plot the CML for dimensions 1 to 3. 
The values of e are chosen as -s/dg x 10% of the fluctuation 
of the time series where d^; is the embedding dimension 
[3^ . The embedding dimension for the CMLs is 3 and 
the graphical analyses are not altered when this number 



is changed. The criterion of nonstationarity is that the 
density of points in the upper left and lower right corners 
of the recurrence plot should decrease ^35^ . As we can see 
from the plots, we obtain inconclusive result in 3D but 
for the lower dimension lattices the system, according to 
this criterion, are stationary. 

The method discussed in Ref. [s^ uses the property 
of nonlinear dynamical systems that the recurrence time, 
suitably calculated, is roughly constant for stationary dy- 
namics. The idea is to compute the recurrence time for 
a trajectory given an initial condition inside a ball of ra- 
dius r > (more specifically, recurrence time of second 
type). We consider an embedding dimension 2 but the 
results are robust with respect to changing this parame- 
ter. In order to eliminate the dependency on the initial 
condition a certain procedure is needed. This requires 
normalization of the recurrent time and then several ini- 
tial conditions are taken. A ball of radius r is centered at 
each initial condition an the times for the first return to 
each ball are computed. The first return time is the aver- 
age of these. Then the second, third, etc, return times are 
computed. The criterion for stationarity is that the plot 
of the time necessary for returns should be constant as 
a function of the number of recurrences. One concludes, 
according to this method, that the CMLs are stationary. 

Another test for stationarity §3] requires the parti- 
tion of the time series and the implementation of some 
predictability algorithm (or any other statistic sensitive 
to differences in dynamics). The algorithm is used as a 
criterion to evaluate how the data from subsets of the 
partition can be used to predict each other (thus the 
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name cross prediction used by Schreiber [3J| ) . Here again 
the choice of embedding dimension does not affect the 
main conclusions and we choose it to be 3, with unit de- 
lay. Suppose the time series Xn, n = 1, 2, . . . , TV, is split 
into contiguous segments of size I and that a statistic 
7y is constructed for each pair of segments For 
the present purpose we choose a simple prediction algo- 
rithm on state vectors [37[ in which 7.^ is represented 
as an error surface on the grid that is, the error 

of using segment i to predict the dynamics in segment 
j. A plot of this surface gives an idea of whether dif- 
ferent dynamics are at operation in different segments: 
when the prediction error is small then segments have 
similar dynamics. For nonstationary systems we would 
expect that the error surface has a systematic low value 
at the diagonal since a segment is being used to predict 
its own dynamics (in-sample prediction). For stationary 
systems the dynamics should be similar in all segments 
and a rough error surface is the result. From several error 
surfaces analyzed there is no evidence of nonstationarity 
in the CMLs. 

Apart from the above, recurrence plots can also be used 
to provide information about the maximal Lyapunov ex- 
ponent. It was recently shown by Letellier ^32i | that recur- 
rence plot analysis provide an estimate of the Shannon 
entropy and thereby indicates the existence of a posi- 
tive maximal Lyapunov exponent. The key idea is to use 
the frequency of occurrence of diagonal segments of non- 
recurrent points in the recurrence plots to estimate the 
Shannon entropy. In this study, we have computed the 
Shannon entropy using the recurrence plots for the CML 
and all the three cases. We find that the estimate turns 
out to be a positive quantity indicating the time series 
obtained in the ID, 2D and 3D CML are chaotic. 



V. SUMMARY AND CONCLUSIONS 

In this paper, we have made a comparative analysis on 
the chaotic times series generated from the one-, two- and 
three-dimensional coupled map lattices based on certain 
commonly employed tools such as fluctuation analysis. 



power spectral density, local dimension and recurrence 
quantification analysis. 

We found that the plot of the root mean square of 
the fluctuation for the 3D lattice shows an anomalous 
behavior so that no scaling can be found for n less than 
about 400. However, there are large scaling regions for 
ID and 2D maps. Further, there is notable increase in the 
persistence if we increase the dimension. In other words 
the three-dimensional CML is more persistent than one- 
and two-dimensional CML. Other important properties 
analyzed reveal that the CMLs are low dimensional and, 
in the 3D case, show a clear indication of 1/f noise. All 
the simulations were replicated a different points of the 
lattices and the conclusions are robust in this respect. 

If we analyze the time series as whole, then the tests 
in Section llV CI show that the system is stationary. How- 
ever, it is evident from the plot of the time series in 
Fig. [He) that there must be some change of the vari- 
ance along the evolution. In order to analyze further 
this question we used the structure function in Section 
nil CI and found that there is stationary behavior only for 
time values larger than about 1000. Thus the nonsata- 
tionarity breaks the scaling for the 3D lattice dynamics 
for sufficiently low values of n. As a next step we in- 
tend to perform extensive studies on the predictability of 
the time series generated by the CMLs using algorithms 
specifically designed for low dimensional chaotic systems 
and also other parameter regimes of the coupled map 
lattices where it exhibits spatiotemporal intermittency. 
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